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ABSTRACT 

The Starlette satellite, launched in February 1975 by the French Centre 
National d' Etudes %>atiales, was designed to minimize the effects of non- 
gravitational forces and to obtain the highest possible accuracy for laser range 
measurements. Analyses of the first four months of ^obal laser tracking data 
have confirmed the stability of the orbit and the precision to which the satellite’s 
position can be established. 

Initial orbit computations using the GSFC GEM-7 gravity model produced 
rms fits of about 8 to 10 meters for arc lengths of 5 days. Throi^L a series of 
gravity model improvements these rms fits have been reduced to 1 to 2 meters 
for the 5 day arcs. An rms fit of 4. 3 meters was obtained for a 90 day arc. 

Five day arcs overlapped by 2. 5 days showed rms satellite position differences 
gei^rally less than 2 meters. Prediction errors at the end of two months were 
less than 30 milliseconds. 
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PRECISION ORBIT COMPUTATIONS FOR STARLETTE 


INTRODUCTION 

On February 6, 1975 the Storiette salellite was launched by the French 
Centre National d' Etudes S^atiales (CNES). Tracking and orbit computation ac- 
curacies since achieved indicate that this satellite is of major significance to 
geophysical and geodetic analyses. 

Starlette (CNES/GRGS, 1975) was desigi^d to obtain the highest possible 
accuracy for laser range measurements and to minimize the effects of non- 
gravitational forces. The satellite was designed to be a high density specularly 
reflecting sphere having a radius of 12 cm and a weight of 47. 295 kg. The low 
area to mass ratio was achieved by using a core which is largely uranium 238. 
The skin is an aluminum alloy containing a total of 60 laser corner reflectors. 
The orbit selected has an inclination of 49. 8° and is well above the high dr 
region, havli^ normal apogee and perigee heights of 1105 and 810 km, 
respectively. 

To determine the contribution Starlette might provide to our studies, w'e 
began an in-depth analysis of all available precision reduced laser data in the 
April to August 1975 period. The data were from the CNES, Smithsonian Astro- 
physical Observatory (SAO), and Goddard ^ace Flight Center (GSFC) systems. 
The objectives of the analyses were to assess the stability and accuracy of the 
orbit, the quality of the laser tracking data, the accuracy of our presently adcpted 
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set of geodetic constants and force models, and the amount of improvement 
which could be readily accomplished by using these data to improve our esti- 
mate of the coefficients in the models. Ihe results of these analyses to date are 
presented in detail in the following sections. 

LASER TRACKING DATA 

A total of 809 passes of precision reduced laser data were available for the 
period of March through August 1975. A large percentage of these data were 
recorded by the SAG stations. In June the total of nine stations tracking con- 
tributed typically eight to ten passes of data per day. On numercHis occasions 
the SAO stations at Arequipa, Peru and Olifantsfontein Republic of S. Africa 
observed three passes per day. Thus, a substantial amount of laser data is 
available. 

The GSFC data have been preprocessed using a program developed by 
Carpenter (1976). Extensive analyses of these data by Carpenter during the 
preprocessing stage have indicated that the precision is on the order of 10 cm 
or better. While we have not conducted an extensive short arc analysis of the 
data from the other stations, the internal consistency of the data in long arc 
orbits from Arequipa, Olifantsfontein, Mt. Hopkins, Arizona and Natal, Brazil 
is on the order of 50 to 150 cm, and is typically 100 to 200 cm for the Greece 
station and the CNES station in the Canary Islands. 
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Data rates for the respective systems are; CNES - every four seconds, 
SAG - every 7. 5 seconds except for the SAO/Natlonal Technical University of 
Athens laser - once per 30 seconds, GSFC - once per second. A typical set of 
residuals for a GSFC site is presented in Figure 1, 

ORBITAL PERTURBATIONS 

To evaluate the nominal behavior of the Starlette orbit, both an analytic 
gravitational harmonic analysis approach and a long periodic semi-analytic ap- 
proach were used. 

The gravitational harmonic analysis approach, embodied as the Harmonic 
Analysis Program (HAP) computer program, computes the expected perturba- 
tion amplitudes associated with each specific frequency in the orbit (i.e. , argu- 
ment of perigee, right ascension of the node, mean anomaly). The formulation 
is taken from Kaula (1966, Eq. 3.76). This is a Unear theory based on analytic 
integration of the Lagrange Planetary Equations over a cycle of th^ particular 
orbit-referred driving frequency under the assumption that the elements and 
rates are constant over the interval, HAP also assumes tl^ J 2 secular rates 
sufficiently represent the total rates. The program uses either a specific grav- 
ity field or Kaula's "rule of thumb" as an estimate of the coefficient values: 

This rule of thumb was develc^d by Kaula from auto-covariance analysis of 
gravimetry and has proved remarkably successful. 


3 


The rates on Storiette indicate that the orbit has a Mrly short perigee 
period of about 110 days, with the observable periods of the zonal perturbations 
being around 110, 55, 37, and 27 days. The nodal period is also fairly short, 
being around 91 days. The nodal period with respect to the sun is about 73 days. 
Taking into account orbit plane symmetry, we would expect to observe solar 
perturbations with a frequency of around 73/2 37 days which is close to the 

iql = 3 zonal frequency and hence with certain kinds of analyses, drag or zonal 
error will resemble solar perturbations. 

These HAP analyses indicate that the primary resonance perturbations have 
beat periods around 2. 7 days due to 14‘*i order terms of the geopotential. When 
Kaula's rule is used to approximate the geopotential, an along track effect of 
about 50 to 100 meters is predicted for these 14**^ order terms. Similar results 
are obtained when the GEM-7 geopotential coefficients (Wagner et al. , 1976) are 
used. In suldition to the 14^^ order resonance, there is a sideband resonance 
with 13^^ order terms having a beat period of around 1. 5 days whose efiect along 
track is in excess of 10 m. The possibility also exists for being able to observe 
the double resonance with 28**^ order terms which have beat periods of around 
1.3 days. 

The non-resonant perturbations are generally less than one meter along 
track for tlie coefficients above (12, 12). However, some higher degree coef- 
ficients, for example (14, 10), produce effects larger than a meter. Even above 
(15, 15), 10cm effects are typically observed. 
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In order to evaluate the pertuxi^ations due to .itmospi^ric drag, direct solar 
radiation pressure, lunar and solar solid body tides and the lunar and solar 
direct gra\'ltational effects, we employed a long periodic semi-analytic program, 
the "Hapid Orbit Analysis and Determination" program {Wagner et al. , 1974). 
ROAD numerically integrates the loz^ periodic orbit and variational equations. 
The ROAD formulation of the potential is similar to that presented in Kaula 
(1966) with die short period variations removed. The above perturbations were 
investigated by comparing orbits with aiui without the desired perturbation o\^r 
the time span for which we had data. The initial conditions are Identical so 
there will be some slight energy discrepancy. 

Figures 2, 3, and 4 present the computed effects of these forces on the 
semi-major axis. Inclination and the right ascension of the ascending node. The 
direct luni-solar perturbations are qualitatively like the tides except for a scale 
factor and thus were intentionally omitted from the graphs. 

Over the three month time period, drag reduces the semi-major axis by 
about 1. 5 meters while radiation pressure causes an increase of about 0. 5 
meters. Radiation pressure produces a maximum perturbation of 0. 02 arc- 
seconds in the inclination while the drag effects of less than 10~^ arcseconds 
are insignificant. % comparison, the tidal perturbations are on the order ;^f 
anarcsecond. In .'he node, radiation pressure and drag produce perturbations 
of about equal magnitude (0. 4 arcseconds) but with 0 {:posite signs. Tidal per- 
turbations in the m>de are on the ortter 1 to 3 arcsecoiKis. 


3 


As shown in Figures 2» 3 and 4, Uie ^rhiiiiing forces produce very difiterent 
effects in the semi-major axis, inclination and node. This fact facilitates sepa- 
rability of the recovered effects when the difiterent elements are used together 
in a solution* Over a longer Ume span the periodic nature of the solar radiation 
pressure perturbation would become evident. 

Errors in the zonal harmonic coefficients may also produce long periodic 
errors in the node. An assessment of the magnitude of this error source is 
required. 

ORBIT COMPUTATIONS 

For detailed analysis of the Starlette laser data, we used the GEODYN 
(Martin, 1972) orbit and geodetic parameter estimation computer program. The 
program employs an order Cowell numerical integration procedure. Satel- 
lite disturbing forces modeled in the program included the gravl^ field of tl^ 
earth in the form of spherical harmonic coefficients, third-body gravitation, 
atmospheric drag, luni-solar solid-earth tides and direct solar radiation pres- 
sure. For the lur ar and solar positions required in the computation of the third- 
body gravitational effects, the Jet Propulsion Laboratory ephemeris DE-69 was 
used. The atmospheric density was modeled using either the Jacchia model at- 
mospheres (Jacchia 1965, 1970, 1971, Jacchia et al, 1968). 

For the initial analyses, the data set was divided into consecutive 5 day arcs. 
This arc length was selected because it corresponds approximately to two resonant 
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beat periods and woiild enable us to investigate the adequacy of the resonant co- 
efficient modeling. The orbit determinations vt;ro performed using die GEM-7 
and GEM-8 gravity models (Wagner, et al. , 1976), the GSFC 1973 station coor- 
dinates (Marsh, et al. , 1973; , aiyd the 1965 Jacchia atmosphere. 

Rms fits as lax'ge as 16 meters vere observed when the GEM-8 model was 
used and almost 11 meters when GEM-7 was used, with both models demonstrat- 
ing a considerable variation from arc to arc. Analysis of die pass by pass tim- 
ing biases for some sample arcs clearly demonstrated the resonance period of 
about 2. 7 days. 

Resonant coefficients for 13^*^, 14^*’ and 28*^ order terms in the GEM-7 
model W 2 re adjusted in a multi-arc solution with 19 cr'nsecutive 5 day arcs in 
the March to June period. The rms fit for each arc was generally reduced by a 
factor of two or more to an overall level of about 4 meters. The fits were also 
significandy more stable from arc to arc. 

For a measure of orbit error other than the regression statistics, w'e com- 
puted radial, cross-track and along-track position differences between the orbits 
associated with each gravity model. The average rms position dificrcnces be- 
tween the GEM-7 and GEM-8 orbits were 3 to 4 meters for radial and cross- 
track components and idxiut 13 meters along-track. A comparison of orbits pro- 
duced by the GEM-7 model containing the adjusted resonance coefficients with 
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those obtained from GEM-7 showed effects of similar magnitude. Tl^ rms dif- 
ferences between the modified GEM-7 orbits and tl^ GEM-8 orbits were larger 
by almost 50-4 

The overall level of Uk: rms of fit produced in the data reductions using the 
GEM-7 modified gravity model was too large to be e?q)lained by anything other 
than geopotentiaJ. mismodelit^. Fortunately, the global distribution of high qual- 
ity laser data well distributed over a full revolution of the argument of perigee 
provided the oj^ortunity to derive a gravity model which was tailored to Starlette. 

combining the normal equatiOwS produced from 26 five day Starlette arcs with 
the complete GEM-7 normal equations, we were able to adjust a complete GEM-V 
size field. The model chosen after several investigations has beer designated 
as PGS-ST4. Table 1 presents rms fits for several five day arcs in June 1975 
which are typical of the level of tit of the tailored model. VVTiereas the GEM-7 
niodel produced rms fits on the order of 10 meters, tlK PGS-ST4 model produced 
rms fits of 1 to 2 meters. 

To assess the relative consistency of the GEM-7 and PGS-ST4 orbits, the 
technique of overlapping orbits was employed. Table 2 presents rms radial, 
cross-track and along-track orbit differences for five day orbits overlapped by 
two and a half days. The rms along-track differences for the PGS-ST4 orbits 
arc less than 2 meters while differences us lai'ge as 18 meters are noted for 
the GEM-7 model. Radially, Uie consistency of tl^ orbits has been considerably 
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improved from a maximum of 3 meters to less than 40 cm. The cross-track 
differences viere also reduced when tte PGS-ST4 model was used but not as 
substantially. 

Because the PGS-ST4 gravity model was derived using five day arcs and 
the previously discussed tests also used five day arcs* arcs in loi^ths of two, 
thirty -five, sixty-six and ninety days were also evaluated. For these investi- 
gations, the Jacchia 1971 mo^l atmosphere was used. Drag and solar radia- 
tion pressure coefficients were adjusted in the thirty-five day and longer arcs. 
The ctegradation in the rms fit as a function of the increased arc lengths was 
approximately linear as is shown in Figure 5, reaching a value of only 4.3 
meters after ninety days. 

The ninety day arc is probably of most interest because it represents a full 
period of the dominant unmodeled perturbation (the K, ocean tide). Table 3 pre- 
sents a breakdown of the statistics by station for this arc. There are a total of 
over 540 passes of data from nine stations. TIk effects of the different system 
noise levels in the SAG and GSFC data is clearly evident in the statistics. To 
demonstrate the consistency in this ninety day orbit, timing errors and range 
biases were calculated from the residuals for each pass of data. Figures 6 and 
7 present the results of these calculations for the Arequipa station which had 
very good coverage throughout the arc. Timing errors were generally less than 
2 milliseconds and the range biases were less than five meters. Plots for the 
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other stations are similar in character. Hie causes for these residual patterns 
are currently under investigation. 

Orbit prediction capability is also of great importance for the laser geo^y- 
namics satellites. The GSFC lasers have narrow beam widths, 1/3 milllradian, 
and acquisition is difficult if the predictions contain large errors. Furthermore, 
in the past it has been necessary to update Uie predictions on a weekly basis for 
satellites such as BE-C, GEOS-1, GEOS-2, and GEOS-3. 

To test the prediction c^abilitLes of the PGS-ST4 gravity model, timing 
error, nd range biases were investigated for a five day orbit propagated over 
a two month period. Initially the orbit was propagated using nominal values for 
drag and solar radiation pressure. At the end of the two month period the tim- 
ii^ error calculated at Arequipa reached a value of nearly 200 milliseconds. 
Using more representative drag and solar radiation pressure values (derived 
from a 35 day arc), the timing errors at the end of the period were reduced to 
about 30 milliseconds and the range biases uere less than 15 meters. Plots of 
these timing errors and the range biases over this period are shown in Figures 
8 and 9. This small prediction error is unparalleled in our previous experience. 
For example, prediction error for the BE-C satellite at the end of 60 days would 
be on the order of seconds. 
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CONCLUSIONS 


These analyses of the Storiette laser tracking data have shown that by sim- 
ply improving the graviiy model, rms fits of a few meters can be achieved for 
arc lengths as long as nii^ty days. The consistency of five day orbits is a few 
tens of centimeters radially and less than two meters along-track. Prediction 
errors at tl^ end of two months are reducible to uiuier 30 milliseconds. Any 
significant improvement will necessarily come from model improvement in fije 
area of tides, polar motion, and other subtle geodynamic effects. 

Continued laser tracking data are expected to provide a strong complement 
to the GEOS- III tracking data for a global center of mass adjustment of the track' 
ing station coordinates. Analyses of the long term evolution of the orbit should 
provide significant contributions to studies of earth and ocean tides. As more 
data are collected, contributions to tectonic motion, polar motion and earth 
rotation studies should be realized. 
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Figure 1. Starlette Laser Range Residuals Single Pass 
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Figure 2. Starlette Semi-Major Axis Perturbations 



Figure 3. Starlette Inclination Perturbations Based Upon Love Number = 0.3 
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Figure 4, Starlette Node Perturbations 


L6 



17 




Figure 6. Starlette Ninety Day Orbit Timing Errors 
At The Arequipa, Peru Laser Station 
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JUNE 26. 1975 


Figure 7, Starlette Ninety Day Orbit Range Biases at the Arequlpa, 

Peru Laser Station 
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ELAPSED TIME (in days) 
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Figure 8. Starlette Orbit Prediction Based Upon A Five Day Arc 
Timing Errors At The Arequipa, Peru Laser Staticm 
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ELAPSED TIME (in days) 

June 26, 1975 

Figure 9. Starlette Orbit Prediction Based Upon a Five Day Arc 
Range Biases At The Arequlpa, Peru Laser Station 
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Table 1 

Starlette Gra’.'ity Model Comparison RMS Fits for 5-Day Orbits 


Epoch 

(1975) 

No, of 

Observations ♦ 

RMS Fit (Meters) 

GEM-7 

PGS-ST4 

June 11 

1500 

9.5 

1.1 

June 16 

2300 

7.6 

1.1 

June 21 

2000 

10,5 

1,7 

June 26 

1200 

10.5 

2.1 


♦AK)roximate, 


Table 2 

Gravity Model Effects on Ortiit Overlap RMS Differences (Meters) 
(5-Day Orbits With a 2. 5 IXiy Overlap) 


Epoch 

(1975) 

Radial 

Cross Track 

Al(mg 

Track 

GEM-7 

PGS-ST4 

GEM-7 

PGS-ST4 

GEM-7 

PGS-ST4 

June 11 

1.4 

0.4 

1.6 

1.5 

9.6 

1.8 

June 16 

1.1 

0.2 

2.5 

0.2 

4.8 

0.9 

June 21 

0. 8 

0.2 

1.6 

0.9 

4.1 

1.4 

June 26 

3. 0 

0.2 

1 2 

0.5 

18.4 

2.0 
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Statistics by Station for Starlette 90 Day Arc 
May 1 to August 1, 1975 
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